The Spectrum of Electromagnetic Jets from Kerr Black 
Holes and Naked Singularities in the Teukolsky 
Perturbation Theory 



Denitsa Staicova ^ and Plamen Fiziev ^'^ 

dstaicova@phys.uni-sofia.bg 
fiziev@phys.uni-sofia.bg 

Abstract 

We give a new theoretical basis for examination of 
the presence of the Kerr black hole (KBH) or the Kerr 
naked singularity (KNS) in the central engine of dif- 
ferent astrophysical objects around which astrophysical 
jets are typically formed: X-ray binary systems, gamma 
ray bursts (GRBs), active galactic nuclei (AGN), etc. 

Our method is based on the study of the exact solu- 
tions of the Teukolsky master equation for electromag- 
netic perturbations of the Kerr metric. By imposing 
original boundary conditions on the solutions so that 
they describe a coUimated electromagnetic outflow, we 
obtain the spectra of possible primary jets of radia- 
tion, introduced here for the first time. The theoretical 
spectra of primary electromagnetic jets are calculated 
numerically. 

Our main result is a detailed description of the qual- 
itative change of the behavior of primary electromag- 
netic jet frequencies under the transition from the KBH 
to the KNS, considered here as a bifurcation of the Kerr 
metric. We show that quite surprisingly the novel spec- 
tra describe linearly stable primary electromagnetic jets 
from both the KBH and the KNS. Numerical investi- 
gation of the dependence of these primary jet spectra 
on the rotation of the Kerr metric is presented and dis- 
cussed. 

Keywords Jets, GRB, Black holes. Naked singulari- 
ties, Teukolsky master equation 
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1 Introduction 

The formation of astrophysical jets is one of the fron- 
tiers in our understanding of physics of compact mas- 
sive objects. CoUimated outfiows are observed in many 
types of objects: brown dwarfs, young stars. X-ray 
binary systems (Mundt et al. (2009)), neutron stars 
(NS), GRBs, AGN (for a brief review see Livio (2004)), 
even galactic clusters, and observational data con- 
tinue to surprise us with the universality of jet for- 
mation in the Universe. For example, according to 
F(^rnii/LAT Collaboration (2009), observations surpris- 
ingly hint for independence of AGN jet-formation from 
the type of their host galaxy. Despite the differences in 
the temporal and dimensional scales and the suggested 
emission mechanisms of the jets, evidences imply that 
there is some kind of common process serving as an en- 
gine for their formation. The true nature of this engine, 
however, remains elusive. 

An example of the difficulties in front of models 
of an engine like this, can be seen in GRB physics. 
GRBs are cosmic explosions whose jets are highly col- 
limated {Oj^t ~ 2° — 5°), highly variable in time (~ 
seconds), extreme in their energy output (~ 10^^ — 
Iq5Q gj.g-| appearing on different redshifts (cur- 
rently observed up to z^^^ = 8.2) ( Mcszaros (2001), 
Zhang & Mcszaros (2002), Burrows (2005), and for a 
review on the online GRB repository see E\-aus et al. 
(2007)). Although there are already several hundred 
GRBs observed by current missions and a lot of details 
are well studied, the theory of GRBs remains essentially 
incomplete, since we lack both an undisputable model 
of their central engine and a mechanism of creation of 
the corresponding jets. 

Currently, GRBs are divided into two classes - short 
and long - based on their duration (with a limit Tgo ~ 
2s, for details see Ghirlanda et al. (2009)). This tem- 
poral division seemed to correlate with their origin and 
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spectral characteristics, leading to the theory that short 
GRBs are likely products of binary mergers of com- 
pact objects (black hole (BH) and neutron star (NS) 
or NS and NS), while long GRBs are an outcome of 
the collapse of massive stars^. This clear distinction, 
however, is questioned by observations of bursts with 
common for the two classes properties (for the lat- 
est observation of such GRB: Antonclli et al. (2009), 
for a review Ripa et al. (2009)) suggesting the exis- 
tence of an intermediate, third class of GRBs and by 
statistical analysis suggesting that the two classes are 
not qualitatively different (Nysewander et al. (2009), 
Shahmoradi & Nemiroff (2009) and Lii et al. (2010)) - 
both evidences of common producing mechanism for 
the two classes. Furthermore, the possibility to explain 
short GRBs by merger of binary systems of BH-NS or 
NS-NS was recently questioned by the lack of detec- 
tion of gravitational waves in 22 cases of short GRBs 
(LIGO Collaboration & Virgo Collaboration (2U1U), 
LIGO Collaboration & Hurley (2007)). At present we 
are not able to explain definitely these observational 
results. Either we have no binary mergers of compact 
objects in sGRB's, or the very large distances to the 
observed events prevent observation of corresponding 
gravitational waves by the existing detectors. The ab- 
sence of gravitational waves from merger of binaries, 
if confirmed by future careful investigations, may seri- 
ously challenge the current paradigm in the sGRB field. 

The dominating theoretical model in the GRB 
physics describing the evolution of already formed jet 
is the "fireball model" introduced by Piran (1999). 
In this model series of shells of relativistic particles are 
accelerated by unknown massive object called the "cen- 
tral engine". The observed lightcurve is produced by 
internal collisions of the shells and by their propaga- 
tion in the circumburst medium. This model, however, 
doesn't explain the process behind the formation of the 
shells or the coUimation of the jet. 

The discrepancies between this model and the obser- 
vational data were outlined recently in Lyutikov (2009), 
most important of them - the yet unexplained afterglow 
decay plateaus (and the sharp decay of the light curve 
afterwards) and the discovered by the mission SWIFT 
flares (additional sharp maxima superimposed over the 
decaying afterglow lightcurve) observed in different 
epochs of the burst and with different fluences (for a 
detailed study of flares see Maxham & Zhang (2009)). 



^The major evidence for the eoUapsar model, the connection of 
long GRBs with supernovae (SN) was confirmed by the obser- 
vations of GRB 980425/SN 1998bw and GRB060218/SN2006aj 
(Mirabal et al, (2006)) - events that started as a GRB and then 
evolved spectrally into an SN. But the question remains why not 
every long GRB is accompanied by SN. 



Although both plateaus and flares imply energy injec- 
tions by the central engine, the flares that may appear 
up to lO^s after the trigger with fluences sometimes 
comparable with that of the prompt emission, question 
the very nature of these energy injections. Different 
theories and simulations are being explored to explain 
these problems (for recent studies see Granot (2006), 
Dado & Dar (2009) , Fan (2009), Shcn et al. (2009), 
Willingale et al. (2009)), but the obvious conclusion 
is that we still lack clear understanding of the central 
engine of GRB. 

Theoretically, the most exploited models of cen- 
tral engines include a Kerr black hole (KBH) in 
super-radiant mode(Whceler (1971), Press & Teukolsky 
(1972)) - the wave analogue of the Penrose process 
or the Blandford-Znajek process (Blandford & Zuajck 
(1977), Blandford (2001)) based on electromagnetic ex- 
traction of energy from the KBH - the electromagnetic 
analogue of Penrose process. Both processes have their 
strong and weak sides, but the main problem is that 
they cannot provide enough energy release in a very 
short time interval, typical of the GRB. The Penrose 
process seems to be not efficient enough (it offers signif- 
icant acceleration only for already relativistic matter), 
as shown long time ago by Wald (1974). This prob- 
lem is also discussed in Bardecn et al. (1972), who ar- 
gue that the energy gap between a bound stable orbit 
around the KBH and an orbit plunging into the ergo 
region is so big, "energy extraction cannot be achieved 
unless hydrodynamical forces or superstrong radiation 
reactions can accelerate fragments to speed more than 
0.5c during the infall". As for the Blandford-Znajek 
process which offers a good explanation of the col- 
limation we observe, it requires extremely intensive 
magnetic fields to accelerates the jets to the energies 
observed in GRBs (theoretical estimations show that 
B - 10^5 10^^ G is needed for a jet with E > 10^^ erg 
) . Even when the magnetic fields are sufficiently strong 
( as in da Silva de Souza & Ophcr (2009) ), numeri- 
cal simulations imply that the BZ process is not effi- 
cient enough (Nagataki (2009)) for the formation of the 
powerful jets seen in GRBs (Punsly & Coroniti (1990), 
Lee et al. (1999), Komissarov (2008), Barkov & Komis- 
sarov: (2008a), (2008b) , Lei et al. (2008), 
Krolik & Hawley (2009), Komissarov & Barkov (2009)). 

Furthermore, the most speculated in the theory GRB 
engines include a rotating black hole, but the observa- 
tional data on GRBs do not provide definitive clarifi- 
cation on this assumption for now. The main problem 
is that we cannot observe the central engine directly 
not only because of the huge distance to the object 
and its small size. The usual techniques to study its 
strong gravitational field are not applicable for the GRB 
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case, because of the properties of their physical pro- 
genitors and their highly non-trivial behavior: 1) We 
are not able to observe the motion of objects in the 
vicinity of the very central engine of GRB. 2) Instead 
of a quiet phase after the hypothetical formation of 
the KBH, we observe late time engine activity (flares) 
which is hardly compatible with the KBH model. 3) 
The visible jets are formed at a distance of 20-100 
event horizon radii (Kocnigl (2006),Vlahakis (2006), 
Marshall ct al. (2006).) where one cannot distinguish 
the exterior field of a Kerr black hole from the exte- 
rior field of another rotating massive compact matter 
object solely by measuring the two parameters of the 
metric - M and a (for a more detailed discussion see 
Fiziev (2009a), Fizicv (2010a)). 

In this situation, the only way to reveal the true na- 
ture of the physical object behind the central engine 
is to study the spectra of perturbations of its gravita- 
tional field. This background field is either described 
exactly by the Kerr metric (Kerr (1963)), if the central 
engine is a KBH or a KNS, or at least approximately 
modeled by the Kerr metric with the same M and a, if 
the central engine is some other object. In any case we 
can use the Kerr metric to study the jet problem in the 
above conditions. 

Different types of central engines yield different 
boundary conditions for the perturbations we intend to 
consider and these boundary conditions generate differ- 
ent spectra specific for the given object. Thus finding 
appropriate spectra to fit our observations enables us 
to uncover the true nature of the central engine. 

Similar method for discovering of the BH horizon 
was proposed for the first time in Dctweilcr (1980), 
and studied more recently in Dreyer et al. (2004), 
Fiziev (2006), Chirenti & Rezzolla (2007), (2008). 

In the present article we give a theoretical basis for 
application of the same ideas for studying the nature 
of the central engine of GRBs, AGN, etc., based on the 
spectra of their jets (Fiziev (2009a)). 

The main new hypotheses studied in the present ar- 
ticle is that the central engine generates powerful pri- 
mary jets of radiation: electromagnetic, gravitational 
and maybe neutrino primary jets. We investigate a 
simple mechanism of formation and collimation of such 
primary jets by the rotating gravitational field of the 
central engine, using a new class of solutions to the 
Teukolsky master equation. Further on, one can imag- 
ine that the primary jets, if powered enough, inject en- 
ergy into the surrounding matter and start the process 
of shock waves and other well studied effects, investi- 
gated in the MHD models of jets, accelerating the ex- 
isting particles, creating pairs of particles-antiparticles 
and so on. This way we are trying to make a new step 



toward filling the existing gap between modern detailed 
MHD numerical simulations of jet dynamics after the 
jets are created and possible mechanisms of the very jet 
creation. The present paper has to be considered as a 
first step in this new direction of investigation. 

Here we examine the electromagnetic-primary-jet 
spectra (see Section 2) of the KBH and the KNS sup- 
posing that the gravitational field of the central engine 
is described exactly by the Kerr metric. 

The study of the spectra of some types of perturba- 
tions of rotating black holes has already serious theo- 
retical and numerical basis, particularly concerning the 
quasi-normal modes (QNM) of black holes. The QNM 
case examines linearized perturbations of the Kerr met- 
ric described by the Teukolsky angular and radial equa- 
tions (the TAE and the TRE accordingly, see in Section 
2) with two independent additional requirements as fol- 
lows: 

A. Boundary condition on the solutions to the TRE: 
only ingoing waves in the horizon, only outgoing waves 
to radial infinity - the so called hlack hole boundary 
condition (BHBC). 

B. Regularity condition is usually imposed on the 
solutions to the TAE (see in Tcukosky (1973) and also 
in Cliaiidrasekliar (1983)). 

The complex frequencies obtained in the case of 
QNM are the solutions of the eigenvalue problem cor- 
responding to these two conditions. They represent the 
"ringing" which governs the behavior of the black hole 
in late epochs and they depend only on the two metric 
parameters M and a. 

A major technical difficulty when searching for QNM 
is that one solves connected problem with two complex 
spectral parameters - the separation constants E and 
oj. This was first done by Teukolsky & Press (1974) 
and later developed through the method of continued 
fractions by Leaver (1985). For more recent results, see 
also Berti ct al. (2004). 

In the recent articles (Fiziev & Staicova (2009a), 
Fizicv & Staicova (2009b)), we illustrate the strength 
of purely gravitational effects in the formation of pri- 
mary jets of radiation. The idea of such physical inter- 
pretation is derived from the visualization of the singu- 
lar solutions to the TAE which shows coUimated out- 
flows. Our toy model is based on linearized electromag- 
netic (spin-one) perturbations of the Kerr black holes 
with modified additional conditions: While we keep the 
BHBC imposed on the TRE, we drop the regularity 
condition on the TAE and we work with specific singu- 
lar solutions of the TAE instead. To achieve this, we 
impose a polynomial condition on the solutions of the 
TAE. This new angular condition reflects the change of 
the physical problem at hand: now we describe a pri- 
mary jet of radiation as a solution to the TAE with an 
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angular singularity on one of the poles of the KBH or 
KNS. 

It is easy to impose the polynomial condition using 
the properties of the confluent Heun functions, which 
enter in the exact analytical solutions of the Teukol- 
sky equations (for more on the use of confluent Heun 
functions in Teukolsky equations see: Fiziev (2009a), 
Fiziev (2009c) and Fiziev (2009d)). 

Here we focus on the case of electromagnetic pertur- 
bations (I s 1= 1 ) which seems to be more relevant 
to the GRB theory, since their observation is due to 
electromagnetic waves of different frequencies - from 
optical to gamma-rays at high energies. Similar situa- 
tion exists for jets from other astrophysical objects. For 
them radio frequencies may be also relevant, say for jets 
from AGNs. In this paper, we elaborate our theoreti- 
cal basis and we report our latest numerical results and 
their theoretical implications. 

Analogous results can be obtained for gravitational 
and neutrino primary jets. Since such jets are still not 
observed in the Nature, we will not include them in the 
present article and we will discuss only electromagnetic 
primary jets. For brevity, from now on we will omit 
"electromagnetic" from electromagnetic primary jets, 
where it doesn't lead to confusion. 

We study in detail the complex spectra of primary- 
jet frequencies and their dependence over the rotational 
parameter of the Kerr metric, a. The parameter a in 
our calculations varies from a = to a = M and then 
from a = M to a >> M. The case a e (0,M) corre- 
sponds to the normally spinning Kerr spacetime; and 
the case a g (M, oo), to the overspinning Kerr space- 
time. 

Much like in flat spacetime, we can consider differ- 
ent physical problems in the Kerr spacetime. Imposing 
BHBC for < a < M we obtain the well studied case 
of the KBH. Under proper boundary conditions which 
have to be specified for a S (M, oo), we obtain different 
physical problems related with the Kerr naked singular- 
ity (KNS). The case a = M under BHBC corresponds 
to the extremal KBH. This case is beyond the scope of 
the present article since it requires a slightly different 
mathematical treatment. 

Our main result is the description of the qualitative 
change of the behavior of the primary jet frequencies 
under the transition from the KBH to the KNS. More- 
over, we discover for the first time that the stability 
condition in our primary jet spectra (positive imagi- 
nary part ensuring damping of the perturbations with 
time) remains fulfilled even in the KNS regime (a > M). 
This happens even though, the imaginary part of those 
frequencies tends to zero as a — ?> M - after this critical 
point, with the increase of a, the imaginary part of the 



frequencies stay positive and grows until it reaches a 
proper constant value for a >> M. 

In our numerical study of primary jet spectra, there 
exist two special frequencies (n = 0, 1) which have real 
parts that coincide with the critical frequency of su- 
perradiance in the QNM case. In our case, the fre- 
quencies, however, are complex, and the imaginary part 
of one of them is of the same order of amplitude as 
the real part, yielding an exponential damping of the 
maybe superradiance-like emission of primary jets, cre- 
ated by the KBH or the KNS. The other frequency 
has a very small (but positive) imaginary part (for 
a < M) . We track the non-trivial change of this special 
superradiance-like frequency with the change of a. The 
relation Wri=o,i,m(") best fitted by formula (3.4a) in 
Fiziev (2009d) for n = TV = 0, 1, m = 0, ±1, ±2 in the 
whole range of a (except for a few points for small a). 
This provides an analytical formula for these modes of 
the spectra. 

As for the other frequencies, they possibly form an 
infinite set as in the case of QNM. We show that the 
primary-jet frequencies differ from the QNM ones due 
to the different condition imposed on the solutions to 
the TAE. The dependence of the primary-jet frequen- 
cies on the parameter a is also presented here for the 
first time. The results are discussed in comparison with 
the case of QNM. 

The article is organized as follows: Section 2 develops 
the mathematical theory behind our spectra. Section 3 
presents the numerical results we obtained, in Section 4 
we summarize our results and we discuss their possible 
physical implications. 

2 A toy mathematical model of formation of 
electromagnetic primary jets by the central 
engine 

The main problem of most GRB models is that they 
do not address the formation of the collimated outflow 
itself, but instead they model the propagation of the jet 
afterwards. 

The propagation of relativistic jets made of parti- 
cles can be best understood by full 3-d general rela- 
tivistic MHD simulations. In numbers of MHD simu- 
lations, however, the energy and/or the angle of col- 
limation of the matter accelerated by the central en- 
gine are just initial parameters (some recent examples 
Mao et al. (2010), Mehani & Keppens (2010)). While 
this approach can bring information about the environ- 
ment of the burst or the properties of the jet, it doesn't 
allow us to study the origin of the jet i.e. the central 
engine. More and more studies try to obtain those pa- 
rameters by evolving real initial conditions like the mass 
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a) 



b) 





Fig. 1 A plot of the ergo surfaces for a) KBH: a < M, b) extremal case: a = M and c) KNS: a > M. One can clearly see 
the ring singularity of the Kerr metric and the change of topology due to the transition from ergospheres to ergo torus 



of the accretion disk/torus around the central object (s) 
or the masses of the binaries in the merger (see for 
example Foucart ct al. (2010), Rezzolla et al. (2010), 
Takcuchi ct al. (2010) ), but there are still many prob- 
lems. The powerful late flares and the huge energies 
of some GRBs, so hard to be fit by the models, force 
us to rethink our understanding of the physics of com- 
pact massive objects. What is clear, however, is that 
precisely such object and its strong gravitational field 
should play a definite role in the central engine of pow- 
erful events like GRBs. Only extreme objects of that 
kind arc able to accelerate and collimatc matter to the 
observed high energies in so short interval of time, using 
the physics we know. 

Trying to uncover possible mechanism able to work 
as a central engine, we wanted to get back to the ba- 
sics - finding the purely gravitational effects that one 
can observe due to a compact massive object. Such 
effects should be intrinsic properties of those objects, 
independent of their environment and they are general 
enough to eventually describe the origin of variety of 
astrophysical jets. 

We started examining the electromagnetic perturba- 
tions of the metric of such object under boundary con- 
ditions different than the already studied in the QNM 
case. Although the full physics of the relativistic jets 
may be outside the scope of the linear perturbation 
theory, understanding the outcome from manipulating 
the boundary conditions in that regime may give us 
important clues about the real picture. Besides, it is 
clear that in the domain of formation of astrophysical 
jets (i.e. 20-100 gravitational radii far from the cen- 
tral engine), the gravitational field is weak enough and 
its small perturbations can be described by linear per- 
turbation theory. Hence, this physics could be at least 
partially unveiled by studying the spectra of radiation 
of that engine in the frame of linear perturbation the- 
ory. 

Thus, in our toy model for formation of primary jets 
of radiation we use a compact massive rotating object - 



a KBH or KNS under appropriate boundary conditions. 
Sufficiently far from the horizon its gravitational field 
is described exactly or at least approximately by the 
Kerr metric. The only quantities on which this metric 
depends are the mass M and the rotational parameter 
a related to the angular momentum J, by a = J /M 
(everywhere in the article c = G = 1). 

For a < M, the Kerr metric has two real hori- 
zons (Kerr (1963)) defined by r^ = M± y/AP-a'^. In 
our studies we will focus on the outer region r > r-|_ 
which is known to be linearly stable with respect to 
QNM. This picture holds for a < M (the normally 
spinning Kerr spacetime). For the overspinning Kerr 
spacctime a > M and r± are complex. When black 
hole boundary conditions are imposed for a < M the 
Kerr metric describes a rotating black hole. The ex- 
tremal limit is reached when a = M and in this case 
r^ — r^ — M. Under proper boundary conditions for 
a > M the Kerr metric describes a rotating KNS. In 
our work, we will discuss only electromagnetic primary 
jets from the KBH or the KNS (s = -1). 

The change of the metric due to the change of the 
parameter a can be seen also in the change of the topol- 
ogy of ergo-surfaces of the Kerr metric (Fig. 1). The 
ergo-surfaces are defined by gtt = 1 — „ii , 3I^'^L'i/a\ = 



which gives the stationary limit surfaces r 

"1: 



E' 



^ where 

= M ± yjAP — cos(6')^. As seen on Fig. 1, for 
a from a < M to a > M we observe a clear transition 
form two ergo-spheres to an ergo-torus. 

From a mathematical point of view, it is more nat- 
ural to work with the dimensionless parameter b = 
M/a £ [0,oo). Obviously the point & = 1 is a bifur- 
cation point for the family of the Kerr ergo-surfaces. 
When the bifurcation parameter b increases from to 
oo, the ergo-torus (6 G (0,1)) bifurcates at the point 
5 = 1 to two ergo-sphercs {b = (1, oo)). The point 6 = 
corresponds to the ring singularity. At the point b = oo 
the outer ergo-sphere transforms to the Schwarzschild 
event horizon and the inner ergo-sphere degenerates to 
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the well-known singularity r = of the Schwarzschild 
metric. 

The evolution of linear perturbations with differ- 
ent spin (I s 1= 0,1/2,1,3/2,2) on the background 
of the Kerr metric was pioneered by Teukolsky in the 
70s and led to the famous Teukolsky Master Equa- 
tion (TME). This differential equation unifies all phys- 
ically interesting perturbations s'^(t,r,9^(j)) written in 
terms of Newman-Pcnrose scalars and is able to de- 
scribe completely multitude of physical phenomena 
(Teukolsky (1972)). 

The TME in the Boyer-Lindquist coordinates is sep- 
arable under the substitution: = e*('^*+™'^)S'(6')i?(r) 
where m = 0,±1,±2... for integer spins and w = 
ujR + iu>i is complex frequency (note that in this ar- 
ticle we use the Chandrasekhar notation in which the 
sign of w is opposite to the one Teukolsky used). This 
frequency and the parameter E are the two complex 
constants of the separation. The stability condition re- 
quires LJi > ensuring that the initial perturbation 
will damp with time. Making the separation, we ob- 
tain the Teukolsky angular (TAE) and the Teukolsky 
radial (TRE) equations which govern the angular and 
the radial evolutions of the perturbation, respectively. 
The key property of those equations is the fact that 
they both can be solved in terms of the confluent Heun 
function as first written in detail in Fiziev (2009a). 

Explicitly, the TAE and the TRE are: 



[auJuY + 2aujsu+ sEir 



{m+suY 



Sim = 0, 



(1) 



and 



+ (1 + .s) 



r — 7 + r — r_ 



dR 



LJ.tj .771 



dr 




(r — r_|-)(r — r_) 



(2) 



a = M , r+ = and the TRE has another irregular 
singularity (r = M) instead of the two regular ones. 
Hence this is a biconfiuent case. Note that the ring 
singularity r = 0, = ■n/2 is not a singularity of the 
TME and does not play any role in its solutions (as 
easily seen from Eq. (1) and Eq. (2)). 

When a 7^ M, these equations can be reduced to 
the confiuent Heun equation and they can be solved in 
terms of confluent Heun's function defined as the unique 
particular solution of following differential equation: 



dz2 



H{z) 



/3 + 1 7+1 



z-l 
z 



d 

dz 

V 



z-l 



H{z) 

H{z) = 0. 



(3) 



The solution called the Heun function is regular in the 
vicinity of the regular singular point z = and is nor- 
malized to be equal to 1 at this point, see the mono- 
graph Slavyanov & Lay (2000) and some additional ref- 
erences in Fiziev (2009a), Fiziev (2009b). 

From the properties of the confluent Heun function 
follows that the general solutions of Eq. (3) can be 
written, using the maple notation, as 

H (z) = Ci HcunC (a, /3, 7, 5, rj, z) + 

C72 2^^HeunC(a,-/?,7,J,7?,z), (4) 

where 5 = ^+i.-a^±2±2 and ry = iii^-^-^±Xtfe. 

According to articles Fiziev (2009a, 2009c, 2009d), 
Fiziev & Staicova (2009a, 2009b) and Borissov & Fiziev 
2009) the solutions of the TAE and the TRE for the 
electromagnetic case s = — 1 can be written in the form: 



[r) = d ''HeunC a, ^, 7, 5, 77, - 



and 



Ri ir) = C2 e''^''HeunC a, -/3, 7, (5, 77 



- + 1 



where A 



2Mr = {r - rj){r - r+), K = 



—uj{r^ + a?) — ma, \ = E — s{s + l) + a^uj^ + 2amuj and 
u = cos{9). 

In general, these equations are ordinary linear 
second-order differential equations with 3 singularities, 
two of which regular, while r = 00 is irregular sin- 
gularity obtained via confluence of two regular singu- 
lar points. For the TRE, the regular singularities r± 
are different when a ^ M. At the bifurcation point. 



in the case, r > r^, where: 

a = —2i (r^ — r_) /? = - 

7 = 1, 



S ~ ~'2i(^r^ — 7'_)a;(l — i (r_ + r^) w), 

(4uj'^r^^+4 (iw-2wV_) r^+{l-4aivm-2uj'^a'^-2E) x 
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4 (iw r_ — 2iLU + E - io"^ — ^) r_ r^—^a^ {m+uo of ^ . 

For the angular function S{9), wc obtain the follow- 
ing solutions: 

HcunC(Q;,/3,7,(5,77,z±(6')) 

where 

a — ±4 aw, 13 = s^m,j = sim, S — =f4 oojs, 

9 9 

77= ±2au}s—(auj) —E+s , 

z+{0) - cos^ (^0 ,z-{e) - sin^ (0=1- 

Note, the confluent Heun function depends on 5 pa- 
rameters (here a, f3,^, 6 and jy ) but the values of those 
parameters differ for the TAB and the TRE. 

To fix the physical problem we want to study, we im- 
pose appropriate boundary conditions on the solutions 
to the TRE and the TAE. These boundary conditions 
yield the spectrum of the separation constants E and 

LU. 

For the solutions of the the TAE, we use the new 
requirement that the confluent Heun functions should 
be polynomial. The polynomiality condition reads: 

- + ^^ + A^+l = 0. Ajv+iA* =0. 

a 2 

It yields a coUimatcd singular solutions Fizicv (2009a). 
Here, the integer iV > is the degree of the polynomial 
and AAr+i(/z) is three-diagonal determinant specified in 
Fizicv (2009a), Fizicv (2009b). 

The polynomial requirement for the angular solu- 
tions fixes the following relation between E and w: 



j,=_i£'^(cj) = —(auj)^ — 2 auj m±2 y (auj)'^ + auj m. (5) 

This simple exact relation leads to a significant tech- 
nical advantage. Instead of working with a complicated 
connected system of spectral equations (as in the QNM 
case), by imposing the polynomial requirement we have 
to solve only one spectral equation for the variable w. 
Thus we essentially simplify the calculations and we ob- 
tain interesting from physical point of view new results. 

Just examining this simple form of the relation Eitj) 
and using the polynomial solutions, we plotted S{6) 
which controls the angular behavior of the solution 
and we were able to obtain different types of col- 
limated outflows, generated by electromagnetic per- 
turbations of the Kerr metric for different values of 
to which not necessarily belong to the jet spectrum 



(Fizicv & Staicova (2009a)). The existence of such col- 
limation in our results justifies the names "jet" solution 
and "jet" spectra that we use. Unlike the QNM case, 
the considered here "jet" case naturally generates a col- 
limated electromagnetic radiation outflow from KBH 
and KNS. In the present article we confirm the above 
preliminary results obtaining numerically the primary- 
jet spectrum of uj defined by the KBH and the KNS 
jet-boundary-conditions. 

It is important to emphasize that in the case of the 
angular equation we work with the singular solutions of 
the differential equation. The question of the possible 
physics behind the use of such singularity is answered 
in Fizicv (2009d). By looking for solution in a specific 
factorized form * = e*("*+""^)S'(6')i?(r), the function 
we obtain in general is not a physical quantity. In- 
stead it defines a factorized kernel sKE.u,m{t, t, 0, </>) ^ 
e'('-t+^n^)^SE^cuM^)sRE,u,rn{r) of the general integral 
representation for the physical solutions of the TME: 

,«-(<, r, 0,0) = 

X] 9^ [di^ [dE sA^,E,7nsKE,uj,ni{t,r,d,(l)). (6) 

This form of the mathematical representation of the 
physical solution is written as the most general super- 
position of all particular solutions of the TME and it as- 
sumes summation over all admissible values of E and uj. 
If one wants to work with the solutions corresponding 
to certain boundary conditions, then one has to account 
for the specific admissible spectra of E and uj in those 
cases. For example, if we have obtained discrete spectra 
for UJ, this leads to the use of singular kernel propor- 
tional to ^„ S{uj — uJn), where 6 is the Dirac function. 
This reduces the integral over uj to summation over aj„ . 
Similarly, for solutions with definite total angular mo- 
mentum we have E = 1(1 + 1) and the integral over E is 
replaced by summation over Z-integer (or half integer), 
because of singular factor J^i — + 1)) Eq. (6). 

From Eq. (6) it becomes clear that the singular so- 
lutions of the angular equation do not cause physical 
difficulties, only if wc are able to find appropriate am- 
plitudes sAu,^E,m that will make the physical solution 
s^{t,r,0,(j)) regular. Examples of finding such ampli- 
tudes are being developed. It is clear that this formula 
justifies the use of singular solutions of the TAE (for 
s = — 1 the singularity is on one of the poles of the 
sphere, for s = 1 - on the other one). 

Knowing E(uj), we can find numerically the frequen- 
cies UJ by imposing boundary condition on the solutions 
to the radial equation. The Kerr spacetime can be con- 
sidered as a background for different physical problems. 
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As in flat spacetimes, one has to fix the pliysics impos- 
ing the corresponding boundary conditions. A specific 
pecuharity of Kerr spacetime is that in general the TRE 
(Eq. (2)) has 3 different singular points: r_|_, r_ and oo. 
To fix the correct physical problem, we have to specify 
the boundary conditions on two of them (the asymp- 
totic behavior of the solution at these singular points). 
This represents the so-called central two-point connec- 
tion problem (as described in Sla\'yaiio\- & La\' (20110)). 
In principle, we can impose boundary conditions on dif- 
ferent pairs of singular points: on and r+; on r_ and 
-\-oo; on r_ and —oo; on r+ and -|-oo; on r+ and — oo; 
and on —00 and -l-oo. Each of these choices will fix a 
specific type of physical problem posed in the complex 
domain of variables. Obviously, these different prob- 
lems may have different physical properties, especially, 
with respect to their linear stability. 

Note that choosing a specific central two-point con- 
nection problem we fix the physical problem indepen- 
dently of the values of bifurcation parameter b. This 
way, we are able to study the bifurcation phenomenon 
for the given physical problem in the whole range of the 
parameter b. Since we want to study primary jets from 
the KBH and the KNS, we start by imposing BHBC on 
the solutions to the TRE at the points r+ and oo for the 
first case (a < M), where the regular singularities are 
real. These boundary conditions are physically well mo- 
tivated (Teukolsky (1972), Chandrasckhar (1983)) and 
the central two-point connection problem for the KBH 
has a clear physical meaning. The same central two- 
point connection problem exists in the overspinning 
case, in spite of the fact that its two regular singulari- 
ties r± are complex. We use the same boundary condi- 
tions for the case of primary jets from the KNS, since 
the physics of the problem is defined by them. Thus, 
we transform the BHBC to naked singularity bound- 
ary conditions using an analytical continuation of the 
central two-point connection problem. 

BHBC can be summarized to: 

1. On the horizon (r — r+), we require only ingoing 
(in the horizon) waves. This specifies which of the two 
solutions of the TRE (Ri or R2) works in each interval 
for the frequency uj. In our case, for any integer m, i?i 
is valid for ^(uj) £ (— mr2+,0), R2 is valid in the com- 
plimentary set 3fi(w) G (—00, +00) \ (— mri-(-,0) where 
f7+ = a/2Mr+. For m = 0, the only valid ingoing 
solution in the whole range 3fi(w) S (—00, -l-oo) is i?2- 

2. At infinity (r — > 00) we allow only outgoing waves. 
Explicitly, at infinity we have a linear combination of 
an ingoing (i?<-) and an outgoing (R-^) wave: 



Since this constant is unknown, we find it indirectly 
from: 



R — C<— i?i 



C^R- 



(7) 



If in this equation we set lim = 0, this will elimi- 

nate the second term in Eq. (7). To do this, we com- 
plexify r and uj and choose the direction in the complex 
plane in which this limit tends to zero most quickly. 
This direction turns out to be arg(r) = 37r/2 — arg{u]), 
i.e., r =1 r I g3/2i7r-jarg(cj) connecting the arguments of 
the complexified r and w (Fizicv (2006)). Having fixed 
that, it is enough to solve the spectral condition in the 
form 



lim 



R 

rZ 



= 



(8) 



where C<_, arc unknown constants. In order to 
have only outgoing waves, we need to have = 0. 



in order to completely specify the primary jet spectra 
sUJ^^ln- For primary jets from the KNS we use the same 
spectral condition in the complex domain (see above). 

Since in this paper we calculate only the case s = — 1, 
for primary jets from the KBH and the KNS, we will 
omit the prefix s = — 1 in front of LOri,m and also the 
index "jet". 



3 Numerical results 

In Fiziev & Staicova (2009a, 2009b), we presented some 
preliminary results for the solutions of the TAE. They 
showed that the polynomial angular solutions describe 
coUimated structures from various types for arbitrary 
w, even if the BHBC (Eq. (8)) is not fulfilled. In this 
article we will focus on the solutions of the radial equa- 
tions, thus fixing the specific spectrum cj„^„i for primary 
jets from the KBH and the KNS from Eq. (8). 

3.1 Numerical methods 

To find the zeros of Eq. (8), we use the software package 
MAPLE which currently is the only one able to work with 
Heun's functions. 

The roots presented here are found with modified by 
the team Miiller algorithm, with precision to 13 digits. 
Although the results for different versions of maple 
may vary in some cases due to improvements in the 
numerical algorithm, key values were checked to have at 
least 6 stable digits, in most cases - more than 8 stable 
digits. The parameters are fixed as follows : s = — 1, 
M=l/2,r = 110. This value for the radial variable r is 
chosen so that it represents the actual numerical infinity 
- the closest point at which the frequencies found by 
our numerical method stop changing significantly for 
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-0,05 

Re(o)) 
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Fig. 2 Complex plot of the frequencies aj„ for a — 0. The 
index n = 0, 1, ... increases witli tlie distance to tiie origin. 
For every n > 0, tliere is a pair of 2 roots ujI{ symmetrical 
witii respect to tlie imaginary axis. In tiie case a = 0, tiie 
frequencies are independent of m 



further increase of r (for fixed a, m, n). In study of the 
dependence on a, the step for the rotational parameter 
is = 0.01 for a < M and becomes adaptive for a > 
M. 

In our work, we have studied in detail cases with 
m = 0,±1,2, covering the cases | m |<| s |, | m | = | s |, 
I m |>| s |. The complete set of our numerical results 
is available by request. 

3.2 Summary of the results 

Graphical representation of our results can be seen in 
figures 2-11. 

In the case a ~ when no rotation of metric is 
present we find a set of pairs of frequencies w^'^ with 
equal imaginary parts and real parts symmetrical with 
respect to the imaginary axis {^{uj}^) = —^{uj^) and 
^{uJn) = 3(a;^)). This set (Fig. 2) is maybe infinite 
if the yet unknown next frequencies follow the same 
line as in Fig. 2. The set is independent of azimuthal 
number m as it should be for a ~ 0. Frequencies are 
numbered with n according to their distance from the 
origin | CiJ„ |, starting with n = for the lowest. 

Although Fig. 2 shows the case without rotation, it 
drastically differs from the results obtained by solving 
the Regge- Wheeler Equation for the Schwarzschild met- 
ric in the QNM case ( Chandrasekhar&Dctweile^l975), 
Chandrasekhar (1983), Leaver (1985), Andersson(1992), 
Ferrari (1996), Ferrari (1997), Kokkotas & Schmidt 
(1999), NoUcrt (1999), Bcrti (2004), Fizicv (2006), 
Fiziev (2007), Ferrari & Gualticri (2008), Bcrti ct al. 
(2009)), because of the change of the conditions for 
the solutions of the TAE. 



Using the so found initial frequencies, we track their 
evolution with the change of the rotational parameter 
a. We work with both frequencies in the pairs for each n 
(cj^'^) and wc consider both signs in front of the square 
root in s^^iEi^{uj) (Eq. (5)), which we denote w^^, 
respectively. On the figures, the number of points we 



0.75- 



^ 0.50-071 



0.25- 
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0.5- 



-1.25 -1 -0.75 -0.50 -0.25 



Re(co ) 

^ 0, r 
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Fig. 3 Complex plot of the two lowest modes ui^^q i m('^) 
for m = 0, 1, 2. With red circles we plot lj^^q i with blue 
crosses - <^^=o i mi the curve we mark some values of a. 
It is clear that for m > 0, i^f^ and cjj"^ coincide almost 
everywhere. For m ^ 0, a < M, the symmetry with respect 
to the imaginary axis breaks 
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present is limited by the abilities of the maple numeri- 
cal procedures evaluating the confluent Heun function. 

The relation w„ ^(a) is nontrivial as can easily be 
seen on Fig. 3 and Fig. 4. 

The results we obtained can be summarized as fol- 
lows: 

• The two lowest modes n = 0, 1 we call special and 
they arc the only ones we could trace to a > M . 

• For the special modes, for a < M we observe 
^{'^tm) ~ ^{•^tm) (for m 7^ and a > 0.1) 
while << '^{oJi^^) (for m ^ Q, a > 0.1, 

< 10"''; in the case m = 0, 9(w^q) de- 
creases more slowly). Since the real parts of the two 
modes are almost equal for a < M and the imag- 
inary part of the n = mode is negligible, in our 
figures in that range, we emphasize only the complex 
frequencies ojf^. 

It's important to note that uJq^-^ (for m > Q, a < 
M) has different behavior from oJq^ - it has very 
small real and imaginary parts (compared to all the 
other modes), and since our numerical precision for 
it is low, we would not discuss that mode. For to < 
the situation reverses and 3fi(a;Q",„) w ^{ujf^), while 
becomes the peculiar one. 

• For TO > the real part of the special frequencies 
5i(a;„=o.i,m) decreases steadily with the increase of a, 
until one of the modes has a minimum at the point 
a = M, the other - at a > M, but close to it - 
see Fig. 4 c), e). In the case to = 0, ^{cuifi) has a 
maximum for a < M and then has a minimum for 
a > M (Fig. 4a)). 

• The imaginary part of the special frequencies 3(wi^„i) 
(Fig. 4 b), d), f)) for all the cases stays positive, it 
decreases for a < M and it tends to zero for a = M. 
For a — M, the solutions of the TRE can no longer 
be expressed in terms of confluent Heun functions. 
In this case, one has to use biconfluent Heun func- 
tion (Fizicv (2009d)). We didn't calculate at this 
point. When a > M, 3(a;„^m) increases again un- 
til it reaches almost constant value at high a. 

• In the range a > M , to ^ we obtain two modes 
with different values of the real and imaginary parts, 
but similar behavior. The real parts of the modes 
both have minima, but at different points. The imag- 
inary part of both modes tends to zero as a — > M, 
but for a >> AI they reach different constants with 
Q{uji „i) — 3(^0,™) ~ 1- These two modes appear for 
TO = ±1,±2,±3 (Fig. 3 and Fig. 4). Wc couldn't 
find other modes for a > M using our numerical 
methods. 

• Modes with n > 1 demonstrate highly non trivial 
behavior and strong dependence of the parameter a, 
even though numerically they cannot be traced to 



a > M and in some cases we are able to fix a very 
limited number of points. 

For TO = 0, those modes persistently demonstrate 
signs of loops, (see Fig. 5 and Fig. 7), which seem to 
disappear for m ^ Q (sec Fig. 9, Fig. 10). 

The real parts of the modes with n > 1 seem 
to form a surface whose physical meaning is yet un- 
known, while their imaginary parts split in two with 
increasing of the rotational parameter (Fig. 8, Fig. 10). 

Note that the BHBC requires frequencies with 

^{^n,m) > —mil+ or di{uJn^,n) < (for TO < 0!) to 

be obtained using the solution R2(j'), while frequen- 
cies with < 5R(ci;„.,„) < —mfl+ to be obtained using 
the solution Ri{r). Since working numerically with 
the solution Ri{r) so far has been challenging, we 
present here only frequencies obtained with the solu- 
tion R2{r). Thus, on Fig. 10, only frequencies whose 
real part is above the n = I mode plotted in magenta 
(or below ^{ui) = 0) satisfy BHBC. This, however, is 
the case only for m 0. For to = 0, according to the 
BHBC the solution R2{r) can be used in the whole 
range 3fJ(w„^o) S (— oo,-foo) and all the frequencies 
wc obtained represent ingoing in the horizon modes. 
• The frequencies uj^m ^^'^ ^nm ^'^'^ symmetric for 
TO = (Fig. 6) for 71 > 0, and they coincide in the 
cases TO ^ 0, n = 1 for a > a„, where a^n tends to 
zero for sufficiently big to (for to = 1, = 0.1, for 
TO = 2, Om = 0.04, for TO = 5, Om = 0.02, for to = 10, 
dm — 0). For n > 1, and lOnm do not coincide. 

(a) (b) 




Fig. 6 A particular case of the symmetry observed in the 
frequencies wj^q for a < M. When m 7^ uJn,m a-nd o;"^ 
coincide for a > am 

• There is symmetry in our spectra, which is confirmed 
for TO = ±0,1,2 for all of the modes with n > 
(except at some points): 

where uj^'"^ correspond to the two frequencies from the 
pairs observed for a = (frequencies with positive 
and negative real part for each to, n on Fig. 2). 
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Fig. 4 Dependence of the real and the imaginary part of two lowest modes (n = 0,1) on a for m = 0,1,2. At the 
bifurcation point b = 1 {a = M = 0.5), j.^) — > implying a critical event at this point. Also J?(lj,'|Lq j^.^) = mO,+ 

for a < M, m ^ 0. On the plots loq^ is marked with the dashed line, cjq^ has been ommited. The solid lines and the 
crosses correspond to ijjf^ - they mostly coincide for m 7^ 



This symmetry allows us to study, for example, 
only the case m > if or when it works numerically 
better. 

From the figures, one can clearly see that, for every 
77i, the relation uji m (a) demonstrates systematic be- 
havior characterized by some kind of critical event at 
the bifurcation point b = M/a = 1, where the damp- 
ing of the perturbation tends to zero (3(a;„,m) — >■ 0). 
Physically, this transition corresponds to the change of 
the topology of the ergo-surface visualized on Fig. 1. 

3.3 Analysis of the results 

The spectra we presented demonstrates a clear transi- 
tion at the bifurcation point a = M, matching the tran- 
sition from a rotating black hole to an extremal black 
hole demanded by the theory of black holes. Since we 
can approach the point a = M from both directions, it 
seems natural to speak also of extremal naked singular- 
ity, for a ^ M but a > M. 

From our figures it is clear that the behavior of the 
spectra when a M from both sides is similar - the 
imaginary part of the critical frequency tends to zero 
at those points and quickly rises to a (different) con- 
stant for a ^ M. We couldn't find any frequencies with 



negative imaginary parts in this case; thus, it seems 
that the primary jets from the KBH and the KNS are 
stable even in the regime of naked singularities. This 
differs from the QNM case, in which QNM from naked 
singularities are unstable (Press & Teukolsky (1973)). 

Another example of unstable spectra can be found in 
Dotti ct al. (2008) where completely different bound- 
ary conditions for the overspinning Kerr space-time are 
considered. There, the authors analyzed the central 
two-point connection problem on the singular interval 
r G (— cx3,oo). As we sec, in our case the primary jets 
are stable. There is no contradiction since we study a 
completely different physical problem. Besides, the dis- 
crete levels of the spectra of the primary jets from the 
KBH and the KNS appear to be smooth in the whole 
interval of the bifurcation parameter b, except for the 
bifurcation point 6=1. 

Another surprise is that numerically 3fi(a;„^m) = 
— 77iil-i_ for a < M for modes n = 0, 1. Here = 
a/2Mr^ is the angular velocity of the event horizon. 

It is curious that we can represent the same relation 
in the form 5ft(a;„^m) = —fnujcr, for n = 0,1, where 
the critical frequency of superradiance for QNM 
Wheeler (1971), Zel'dovich (1971), Zel'dovich (1972), 
Starobinskiy (1973), Starobinskiy & Churilov (1973), 
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Fig. 5 Complex plot of some of the frequencies with different n obtained in the case m = 0, a) cj^q, n = 2 — 7, b) cj^ o, 
n = 2 — 7, c) oj^Q for n = 2 — 18 (the black dot stays for a = 0, green diamond - the last a where a root was found). 
The modes n = 0, 1 are not plotted. For m = 0, there is interesting symmetry between figs, a) and b) with respect to the 
vertical axis (see the text for details) 
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a) b) c) d) 

Fig. 7 Two cases showing the detailed behavior of modes with n> 1 for m = Q. a), b) correspond to n = 2; c), d) to n = Q 



Tcukosky (1973), Wald (1974), Kodama (2008) , 
Cardoso ct al. (2004)). 

The complexity of the frequency wi,™ whose real 
part equals —mujcr may show another essential phys- 
ical difference between our primary jets-from-KBH and 
jets-from-KNS solutions and standard QNM, obtained 
using regular solutions of the angular equation Eq. (1). 

In our case, we observe two (lowest) modes {oJq^, 
^fin m > 0) whose real parts coincide with the crit- 
ical frequency of QNM-superradiance. However, only 
'^Om almost real, having a very small imaginary part. 
The magnitude of the imaginary parts of the other two, 
i^fm, is of the same order as that of their real parts. 
To the best of our knowledge this is the first time that 
an imaginary part of such quantity, somehow related 
with the critical QNM-superradiance frequency, is dis- 
covered solving BH boundary conditions in pure vac- 
uum (i.e. without any mirrors, additional fields, etc). 
Note that while for a < M the small imaginary part of 
'^Om speaks of slow damping, for a > M the imaginary 
part seems to grow quickly to become comparable with 
the real part value, thus ensuring stability. The two 



lowest modes in our primary jet spectra arc the only 
ones, which we can trace from a < AI to a > M. Al- 
though the lack of the higher modes (n > 1) for a > M 
could be due to problems with the numerical routines 
we use, its persistence for every m implies that there 
may be a possible deeper physical meaning behind it. 

It is interesting to investigate in detail the relation 
of our results with the superradiance phenomenon. 

Although the considered primary jet-perturbations 
to the Kerr metric in general damp quickly for b =/= 
1 (a 7^ M), this situation changes in the limit b -> 
1 ~- the imaginary part of lo then tends to zero near 
the bifurcation point, both for 6^1 — and b — > 
1 + 0. In this case, the perturbation will damp very 
slowly with time while oscillating violently, opening the 
possibility for interesting new primary jet-phenomena 
which deserve additional detailed consideration. 

The best fit for our numerical data for the lowest 
modes (n = 0, 1) turns out to be the exact formula 
(3.4a) in Fiziev (2009d): 



w„=o,i,m = i-m + iN^b^ N = 0,l. (9) 
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a) h) c) d) 

Fig. 8 The case m — Q, modes with n > 1 a), b): K(cj^o) s^nd Si(ctj^Q), n = 2 — 5 c), d): K(cj^o) and Si(ctj^Q), n = 1 — 18. 
On c) and d) we plotted the conjugate roots a;^'^ on one plot. The magenta lines correspond to the special mode n = 1 




Fig. 9 Complex plots in the case m = —1: a) oj^ n = 3 — 14 , b) ij^ _j, n = 2 — 16, c) ctj,^ n = 2 — 18 . There are no 
loops for j m 1= 1 



Here, the formula is written in terms of angular ve- 
locity r2+ and bifurcation parameter b = M/a since 
these are more relevant to the problem at hand. The 
fit can be seen on Fig. 11. 

Clearly this formula offers a very good fit to the data, 
except for the case to = 0, b > 1 where it gives zeros, 
differing from our numerical results. If we compare all 
other relevant points produced in our numerical evalua- 
tions (i.e., the points for the two lowest n = 0, 1 modes 
for 6 > 1 corresponding to iV = 0, 1 in Eq. (9) and 
points for the only two modes available for 6 < 1, which 
corresponds to = 0, 1 in Eq. (9), for to = 0, ±1, ±2) 
with the values calculated according to formula Eq. (9) 
for the corresponding a and M, the two sets coincide 
with precision of at least 5 digits in the worst few cases, 
and with at least 10 digits for most of the points. Usu- 
ally, there are some deviations for small a (b >> 1). 
The case n = for a < M is not plotted, because nu- 
merically the real parts of n = and n = I coincide for 
0.1 < a < M , while the imaginary part of rt = is very 
small - thus, the n = mode clearly obeys Eq. (9) for 
iV = 0. Taking into account that: 1) This exact for- 
mula is analytically obtained from the properties of the 



confluent Hcun functions in the case of polynomial con- 
dition imposed on the the TRE and not on the TAE; 
and 2) The routines with which MAPLE evaluates the 
Heun function, as well as their precision are unknown, 
this match is extremely encouraging. Besides, it poses 
once more the question about the physics behind Eq. 
(9). 

Eq. (9) appears without derivation also in 
Berti ct al. (2003) and Berti k Kokkotas (2003), where 
the authors imposed QNM boundary conditions mean- 
ing they worked with the regular solutions of the TAE. 
In their paper, Eq. (9) was considered as a linear ap- 
proximation to some unknown nonlinear with respect 
to f2+ relation. It is written in terms of the tempera- 
ture of the horizon, and according to the authors, it fits 
very well the frequencies they obtained for the QNM 
case / = TO = 2 (=| s |) and also it matches the imagi- 
nary parts of the QNM {n = 1, 2, 3...) when to > (for 
the other cases it does not work). 

In contrast, in our numerical results, Eq. (9) de- 
scribes the lowest modes (n = 0, 1) for all values of to 
equally well (except to ~ 0,a < M), while it fails to 
describe the modes with n > 1 - Eq.(9), for iV = n 
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Fig. 11 Comparison between our numerical results plotted with blue (red) crosses and the analytical formula (violet lines) 
in the cases = 0, 1. The excellent fit is spoiled only for m — 0, a < 0.5 



or iV = n + 1 matches only the relative magnitude of 
the modes and has no resemblance to their highly non- 
trivial behavior. The fit does not improve even for the 
highest mode (n = 20) we were able to obtain. Sur- 
prising similarity between the two cases (primary jets 
and QNM) is that the real parts of the higher modes 
(n > 1) in our results also seem to tend to the lowest 
mode 5R(wo,m) = 3?(a;i_m) = — 77iO_|_ . However, since in 
our case the frequencies come in pairs, the frequencies 
with opposite to fi+ sign do not tend to fl^ but they 
still seem to have a definite limit. 

Despite the differences, the appearance of Eq. (9) in 
the two cases corresponding to different physical prob- 
lems is interesting and raises questions about the un- 
derlying physics it represents. It also contributes to the 
validity of our approach, since our results are compa- 
rable in some points with Berti et al. (2003) who use a 
completely different but well established in QNM cal- 



culations numerical approach (Leaver's method of con- 
tinued fractions). 

Finally, it is worth recalling that only part of the 
modes with n > 1, m 7^ presented here satisfy BHBC 
and represent perturbations ingoing into the outer hori- 
zon ?■-(_, thus describing a perturbation of the KBH. 
The modes that violate BHBC (on Fig. 10 between 
3fi(w) = and 3?(wi,_i) - the magenta line) describe 
perturbations outgoing from the horizon, thus describ- 
ing a perturbation of white hole. It is clearly interesting 
that the solutions of the equation R2{r) evolve with 
the rotation a in a way that they first describe ingo- 
ing waves and later describe outgoing waves, with some 
hints that they will tend to fl^, but it is impossible 
to make a conclusion based on current numerical algo- 
rithms and without further study of the whole white 
hole case. Although this is beyond the scope of current 
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article, it clearly demonstrates the complexity of our 
approach and the richness of the possible results. 

4 Conclusion 

In this article, we presented new results of the numer- 
ical studies of our model of primary jets of radiation 
from the central engine inspired by the properties of 
GRB. We have already showed that the polynomial re- 
quirement imposed on the solutions to the TAE leads to 
coUimated jet-like shapes observed in the angular part 
of the solution, thus clearly offering a natural mecha- 
nism for collimation due to purely gravitational effects 
(Fizicv (2010a)). 

Continuing with the TRE, we imposed standard 
black hole boundary conditions and we obtained highly 
nontrivial spectra w„_m(a) for a coUimated electromag- 
netic outflow from both the KBH and the KNS. These 
novel spectra give a new theoretical basis for examina- 
tion of the presence of the Kerr black hole or the Kerr 
naked singularity in the central engine of different ob- 
served astrophysical objects. 

We gave a detailed description of the qualitative 
change of the behavior of primary jet spectra under 
the transition from the KBH to the KNS and showed 
that this transition can be considered as bifurcation of 
the family of ergo-surfaces of the Kerr metric. 

Numerical investigation of the dependence of these 
primary jet spectra on the rotation of the Kerr met- 
ric is presented and discussed for the first time. It is 
based on the study of the corresponding novel primary 
jet-boundary-problem formulated in the present article. 
Here we use the exact solutions of the Teukolsky Mas- 
ter Equation for electromagnetic perturbations of the 
Kerr metric in terms of confluent Heun's function. 

The numerical spectra for the two lowest modes 
n = 0, 1 (with some exceptions) is best described by 
the exact analytical formula in Fizicv (2009d), though 
derived by imposing the polynomial conditions on the 
solutions to the TRE instead to the TAE. Although 
unexpected, this fit serves as confirmation of both our 
numerical approach and the analytical result, and it is a 
hint for deeper physics at work in those cases. Interest- 
ingly, the formula in question, Eq. (9), works only for 
the lowest modes {n = 0, 1) for each m and correctly 
describes two frequencies with coinciding for a < M 
real parts and different imaginary parts (zero and com- 
parable to the real part, respectively). The modes with 
higher n are not described by it and have highly non- 
trivial behavior - they seem to have a definite limit 
which numerically coincides with the angular velocity 
of the event horizon (and critical frequency of QNM 
superradiance) in some of the cases. 



The obtained primary jet spectra seem to describe 
linearly stable electromagnetic primary jets from both 
the KBH and the KNS. Indeed, the imaginary part of 
the complex primary jet frequencies in our numerical 
calculations remains positive, ensuring stability of the 
solutions in direction of time-future infinity and indi- 
cating an explosion in the direction of time-past infin- 
ity. This is surprising, since according to the theory of 
QNM of the KNS, the naked-singularity regime should 
be unstable in time future. There is no contradiction. 
The different conditions on the solutions to the TAE 
wc use, define a different physical situation - we work 
with primary jets from the KBH or the KNS, not with 
QNM, and our numerical results suggest that these pri- 
mary jets are stable even for the KNS. 

We showed that the value b = M/a = 1 describes 
bifurcation of ergo-surfaces of the Kerr metric. Around 
this bifurcation point the imaginary part of the primary 
jet frequency decreases to zero suggesting slower damp- 
ing with time. This could have interesting implications 
(including the loss of primary jet stability) for rotating 
black holes or naked singularities close to the extremal 
regime. Besides, for a >> M (i.e. 6 — >■ 0), or a << M 
(i.e. b — >■ oo) the imaginary part of the primary jet fre- 
quencies remains approximately constant with respect 
to the bifurcation parameter b, being positive in both 
directions. 

Although very simple, our model was able to pro- 
duce interesting results whose physical interpretation 
and applications require additional consideration. It 
also demonstrates the advantage of the use of confiu- 
ent Heun functions and their properties in the physical 
problem at hand. 

Whether this mechanism of collimation is working in 
GRBs, AGNs, quazars and other objects can be under- 
stood only trough comparing the spectra of their real 
jets and our "primary jet spectra" taking into account 
also the interference from the jet environment and from 
the different emission mechanisms. For the real obser- 
vation of the predicted here spectra is important the 
intensity of their lines. It depends on the power of the 
primary jet, created by the central engine, which is still 
an open problem. If there exist a natural gravitational 
mechanism of creation of very powerful primary jets, 
that jets may play the main role in observed astrophys- 
ical jets. Otherwise one can expect to see some traces 
of these primary jet spectra in the spectra of astrophys- 
ical jets, similar to the traces of standard QNM spectra 
in the exact numerical solutions of Einstein equations. 

It is remarkable that similar effects of specific col- 
limation of flows of matter particles, described by 
geodesies in rotating metrics like Kerr one are discov- 
ered independently of our study of specific jet-like solu- 
tions of Teukolsky Master Equation. See the recent 
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articles: Chiconc & Mashhoon (2010), Chiconc ct al. 
(2005) and the references therein. The proper com- 
bination of these novel effects of collimation for both 
radiation and matter outflows from compact objects 
may be important ingredient of the future theory of 
real astrophysical jets. 

The approach based on the Tcukolsky Master Equa- 
tion seems to give more general picture of the collima- 
tion phenomena, at least in the framework of pertur- 
bation theory of wave fields of different spin in Kerr 
metric. In addition, our approach gives direct sugges- 
tions for observations, since the frequencies, we have 
discussed in the present paper, in principle may be 
found in spectra of the real objects like jets from GRBs, 
AGNs, quazars, e.t.c. Their observation will present 
indisputable direct evidences for the real nature of the 
central engine, since the corresponding spectra can be 
considered as fingerprints of that still mysterious ob- 
ject. 
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